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Abstract. The latest development in determining the density dependence of 
the nuclear symmetry energy using heavy-ion collisions is reviewed. Within 
the IBUU04 version of an isospin- and momentum-dependent transport model 
using a modified Gogny effective interaction, recent experimental data from 
NSCL/MSU on isospin diffusion are found to be consistent with a nuclear sym- 
metry energy of E sym (p) « 31.6(p/po) 105 at subnormal densities. Predictions 
on several observables sensitive to the density dependence of the symmetry en- 
ergy at supranormal densities accessible at GSI and the planned Rare Isotope 
Accelerator (RIA) are also made. 
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1. Introduction 

The Equation of State (EOS) of isospin asymmetric nuclear matter can be writ- 
ten within the well-known parabolic approximation, which has been verified by all 
many-body theories, as 



E(p, S) = E{p, 5 = 0)+ E sym (p)S 2 + 0(<5 4 ), 



(1) 
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where S = (p n — p P )/(p P + pn) is the isospin asymmetry and E sym (p) is the density- 
dependent nuclear symmetry energy. The latter is very important for many in- 
teresting astrophysical problems^, the structure of radioactive nuclei E] and 
heavy- ion reactions [0J |SJ |7j. Unfortunately, the density dependence of symmetry 
energy E sym (p), especially at supranormal densities, is still poorly known. Predic- 
tions based on various many-body theories diverge widely at both low and high 
densities. In fact, even the sign of the symmetry energy above 3po remains uncer- 
tain [SJ- Fortunately, heavy-ion reactions, especially those induced by radioactive 
beams, provide a unique opportunity to pin down the density dependence of nuclear 
symmetry energy in terrestrial laboratories. Significant progress in determining the 
symmetry energy at subnormal densities has been made recently both experimen- 
tally and theoretically [ El HHj. High energy radioactive beams to be available at 
GSI and RIA will allow us to determine the symmetry energy at supranormal den- 
sities. In this talk, we highlight the recent most exciting progress in determining 
the symmetry energy at subnormal densities and present our predictions on several 
most sensitive probes of the symmetry energy at supranormal densities. 



2. An isospin- and momentum-dependent transport model for 
nuclear reactions induced by radioactive beams 

Crucial to the extraction of critical information about the E sym (p) is to compare 
experimental data with transport model calculations. We outline here the major 
ingredients of the version IBUU04 of an isospin- and momentum-dependent trans- 
port model for nuclear reactions induced by radioactive beams [ ITT]. The single 
nucleon potential is one of the most important inputs to all transport models. Both 
the isovector (symmetry potential) and isoscalar parts of this potential should be 
momentum dependent due to the non-locality of strong interactions and the Pauli 
exchange effects in many-fermion systems. In the IBUU04, we use a single nucleon 
potential derived from the Hartree-Fock approximation using a modified Gogny 
effective interaction (MDI)[ 12-, i.e., 

U(p,S,p,T,x) = A u (x)— + Ai{x) — 

Po Po 

+ B(-?-y(l-x5 2 )-8TX f/V «W 
Po v + 1 Po 

2C t ,t f 3 , frViP) 



drp 



Po J ^ i + (p-p0 2 /a 2 

f ^ P ' £ (f ff> ■ (2) 



Po J y l + ip-ff/A 2 ' 

In the above r = 1/2 ( — 1/2) for neutrons (protons) and r ^ r'; a — 4/3; f T {r,p) 
is the phase space distribution function at coordinate r and momentum p. The pa- 
rameters A u (x) , Ai (x) , B, C t ,t, C T , T > and A were obtained by fitting the momentum- 
dependence of the U(p,S,p,T,x) to that predicted by the Gogny Hartree-Fock 
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and/or the Brueckner-Hartree-Fock calculations, the saturation properties of sym- 
metric nuclear matter and the symmetry energy of 30 MeV at normal nuclear matter 
density po = 0.16 fm~ 3 [ 2j. The incompressibility Kq of symmetric nuclear matter 
at po is set to be 211 MeV. The parameters A u (x) and Ai(x) depend on the x param- 
eter according to A u (x) = -95.98-a;^- and A L (x) = -120.57+x^-. The param- 
eter x can be adjusted to mimic predictions on the E sym (p) by microscopic and/or 
phenomcnological many-body theories. The last two terms contain the momentum- 
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Fig. 1. Density dependence of the symmetry energy for four x parameters. 



dependence of the single-particle potential. The momentum dependence of the 
symmetry potential stems from the different interaction strength parameters C r ,r' 
and CV )T for a nucleon of isospin r interacting, respectively, with unlike and like 
nucleons in the background fields. More specifically, we use C un iike = —103.4 MeV 
and Cuk e — —11.7 MeV. As an example, shown in Fig. 1 is the density dependence 
of the symmetry energy for x = —2, —1, and 1. 

Systematic analyses of a large number of nucleon-nucleus and (p,n) charge 
exchange scattering experiments at beam energies below about 100 MeV indi- 
cate undoubtedly that the symmetry potential at po, i.e., the Lane potential, 
decreases approximately linearly with increasing beam energy Ekin according to 
U L ane = a- bE km where a ~ 22 - 34 MeV and b ~ 0.1 - 0.2 [EH E]. This provides 
a stringent constraint on the symmetry potential. The potential in eq[3 meets this 
requirement very well as seen in Fig. 2 where the symmetry potential (U n — U p )/25 
as a function of momentum and density for the parameter x = — 1 is displayed. 

One characteristic feature of the momentum dependence of the symmetry poten- 
tial is the different effective masses for neutrons and protons in isospin asymmetric 
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Fig. 2. Symmetry potential as a function of momentum and density for MDI interaction 
with x = — 1. 

nuclear matter, i.e., 



where is the nucleoli Fermi wave number. With the potential in eq. |2 we found 
that the neutron effective mass is higher than the proton effective mass and the 
splitting between them increases with both the density and isospin asymmetry of 
the medium[lll|. 

Since both the incoming current in the initial state and the level density of 
the final state in nucleon-nucleon (NN) scatterings depend on the effective masses 
of colliding nucleons in medium, the in-medium nucleon-nucleon cross sections are 
expected to be reduced by a factor a™]f mm / 'a n n = (^nn/^nn) 2 , where jUaw 
and [i* NN are the reduced mass of the colliding nucleon pairs in free-space and 
in medium, respectively [ E|. Thus, because of the reduced in-medium nucleon 
effective masses and their dependence on the density and isospin asymmetry of the 
medium, the in-medium NN cross sections are not only reduced compared to their 
free-space values, the nn and pp cross sections are split and the difference between 
them grows in more asymmetric matter. We expect the isospin-dependence of the 
in-medium NN cross sections to play an important role in nuclear reactions induced 
by neutron-rich nuclei [H5|. 
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3. Probing the symmetry energy at subnormal densities with 
isospin diffusion 

Tsang et al.[E] recently studied the degree of isospin diffusion in the reaction 124 Sn 
+ 112 Sn by measuring [fTTj 
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(4) 



where X is the average isospin asymmetry (S) of the 124 Sn-like residue. The data 
is indicated in Fig. 3 together with the time evolutions of Ri and the average 
central densities calculated with x = — 1 using both the MDI and the soft Bertsch- 
Das Gupta-Kruse (SBKD) interactions are also shown. It is seen that the isospin 
diffusion process occurs mainly from about 30 fm/c to 80 fm/c corresponding to 
an average central density from about 1.2/9 to 0.3po- The experimental data from 
MSU are seen to be reproduced nicely by the MDI interaction with x = — 1, while 
the SBKD interaction with x — — 1 leads to a significantly lower Ri value [ITU]. 
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Fig. 3. The degree of isospin diffusion as a function of time with the MDI and SBKD 
interactions. The corresponding evolutions of central density are also shown. 

Effects of the symmetry energy on isospin diffusion were also studied by vary- 
ing the parameter x[ I1U| . Only with the parameter x = — 1 the data can be 
well reproduced. The corresponding symmetry energy can be parameterized as 
E S ym{p) ~ 31.6(/9//3o) 105 . In the present study on isospin diffusion, only the free- 
space NN cross sections are used and thus effects completely due to the different 
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density dependence of symmetry energy are investigated. As the next step we are 
currently investigating effects of the in-medium NN cross sections on the isospin 
diffusion. 

4. Probing the symmetry energy at supranormal densities 

Several probes that are sensitive to the high density behavior of the symmetry 
energy have been proposed. As an illustration, we present here three most sensitive 
observables that can be measured in future experiments at R1A and GSI. 

4-1. Pions yields and n~ /n + ratio at RIA and GSI 

At the highest beam energy at RIA, pion production is significant. Pions may thus 
carry interesting information about the EOS of dense neutron-rich matter [ IT5I I19| . 
Shown in Fig. 4 are the 7r~ and 7r + yields as a function of the x parameter. It is 
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Fig. 4. The ir~ and ir + yields as functions of the x parameter. 

interesting to see that the tt~ multiplicity depends more sensitively on the symme- 
try energy. The tt~ multiplicity increases by about 20% while the 7r + multiplicity 
remains about the same when the x parameter is changed from -2 to 1. The mul- 
tiplicity of 7r~ is about 2 to 3 times that of tt + . This is because the ir~ mesons 
are mostly produced from neutron-neutron collisions. Moreover, with the softer 
symmetry energy the high density region is more neutron-rich due to isospin frac- 
tionation[H9j. The ir~ mesons are thus more sensitive to the isospin asymmetry of 
the reaction system and the symmetry energy. However, one should notice that it is 
well known that pion yields are also sensitive to the symmetric part of the nuclear 
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EOS. It is thus hard to get reliable information about the symmetry energy from 
tt~ yields alone. Fortunately, the n~ /ir + ratio is a better probe since statistically 
this ratio is only sensitive to the difference in the chemical potentials for neutrons 
and protons [I20|. This expectation is well demonstrated in Fig. 5. It is seen that 
the pion ratio is quite sensitive to the symmetry energy, especially at low transverse 
momenta. Thus, it is promising that the high density behavior of nuclear symmetry 
energy E sym {p) can be probed using the 7r~ /tt + ratio. 




4-2. Isospin fractionation and n-p differential flow at RIA and GSI 

The degree of isospin equilibration or translucency can be measured by the rapidity 
distribution of nucleon isospin asymmetry Sf ree = (N n — N p )/(N n + N p ) where 
N n (N p ) is the multiplicity of free neutrons (protons) [El . Although it might be 
difficult to measure directly Sf ree because it requires the detection of neutrons, 
similar information can be extracted from ratios of light clusters, such as, 3 H/ 3 He, 
as demonstrated recently within a coalescence model[ l22l 123) . Shown in Fig. 6 are 
the rapidity distributions of 5f ree with (upper window) and without (lower window) 
the Coulomb potential. It is interesting to see that the 5f ree at midrapidity is 
particularly sensitive to the symmetry energy. As the parameter x increases from 
—2 to 1 the 5f ree at midrapidity decreases by about a factor of 2. Moreover, 
the forward-backward asymmetric rapidity distributions of 5/ ree with all four x 
parameters indicates the apparent nuclear translucency during the reaction[ l24j . 
Another observable that is sensitive to the high density behavior of symmetry 
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Fig. 6. Isospin asymmetry of free nucleons with and without the Coulomb force. 

energy is the neutron-proton differential flow[!25j 

N(y) 

i=l 

where Wi — 1(— 1) for neutrons (protons) and N(y) is the total number of free 
nucleons at rapidity y. The differential flow combines constructively effects of the 
symmetry potential on the isospin fractionation and the collective flow. It has the 
advantage of maximizing the effects of the symmetry potential while minimizing 
those of the isoscalar potential. Shown in Fig. 7 is the n-p differential flow for the 
reaction of 132 SYt + 124 Sn at a beam energy of 400 MeV/nucleon and an impact 
parameter of 5 fm. Effects of the symmetry energy are clearly revealed by changing 
the parameter x. 

5. Conclusions 

The density dependence of the symmetry energy is very important for both nu- 
clear physics and astrophysics. Significant progress has been made recently by the 
heavy-ion community in determining the density dependence of the nuclear sym- 
metry energy. Based on transport model calculations, a number of sensitive probes 
of the symmetry energy have been found. The momentum dependence in both the 
isoscalar and isovector parts of the nuclear potential was found to play an impor- 
tant role in extracting accurately the density dependence of the symmetry energy. 
Comparing with recent experimental data on isospin diffusion from NSCL/MSU, 
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Fig. 7. Neutron-proton differential flow at RIA and GSI energies 



we have extracted a symmetry energy of E sym (p) w 31.6(p/ / oo) 1 ' 05 at subnormal 
densities. It would be interesting to compare this conclusion with those extracted 
from studying other observables. More experimental data including neutrons with 
neutron-rich beams in a broad energy range are needed. Looking forward to exper- 
iments at RIA and GSI with high energy radioactive beams, we hope to pin down 
the symmetry energy at supranormal densities in the near future. 
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